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Abstract 

We consider the macroscopic model derived by Degond and Motsch from 
a time-continuous version of the Vicsek model, describing the interaction ori- 
entation in a large number of self-propelled particles. In this article, we study 
the influence of a slight modification at the individual level, letting the re- 
laxation parameter depend on the local density and taking in account some 
anisotropy in the observation kernel (which can model an angle of vision). 

The main result is a certain robustness of this macroscopic limit and of 
the methodology used to derive it. With some adaptations to the concept 
of generalized coUisional invariants, we are able to derive the same system of 
partial differential equations, the only difference being in the definition of the 
coefficients, which depend on the density. This new feature may lead to the 
loss of hyperbolicity in some regimes. 

We provide then a general method which enables us to get asymptotic 
expansions of these coefficients. These expansions shows, in some effective 
situations, that the system is not hyperbolic. This asymptotic study is also 
useful to measure the influence of the angle of vision in the final macroscopic 
model, when the noise is small. 
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1 Introduction 



The study of complex particle systems has given rise to some challenging issues [6], 
in a mathematical point of view. One of the interesting problem is to understand 
how a collective behavior can emerge with only localized interactions. 

The Vicsek model [40] has been proposed as a minimalist model describing the 
behavior of individuals inside animal societies such as fish schools or flocks of birds. 
It is a minimal version of a more complete and realistic model [3, 36, 27, 12] based on 
three zones (of repulsion, alignment, and attraction). The Vicsek model only con- 
siders the alignment behavior, getting around the problem of confining the particles 
in the same region by imposing spatial periodicity (particles move on the fiat torus). 
All the particles have constant speed and synchronously update their direction ac- 
cording to their neighbors, their new orientation vector being given by the mean 
direction (subjected to some angular noise) of all particles at distance less than a 
given radius. As the noise decreases (or the density increases), one can observe a 
phenomenon of phase transition, from a regime of disordered particles, to an ordered 
phase with strong correlations between orientations of particles [40, 2, 24]. 

Two main difficulties arise when we try to derive a macroscopic limit of this 
individual based model. First of all, the system is discrete in time: the time step is 
fixed, and the model is not built in the goal of letting it tend to zero. The second 
problem is that, except for the total mass, there is no obvious conservation relation, 
so a good candidate for a macroscopic model would probably be a non-conservative 
system of partial differential equations, but we lack conservation relations to obtain 
any equation other than the conservation of mass. 

In [15], Degond and Motsch have proposed an approach to handle these two 
complications. First, they provide a time-continuous version of the individual based 
model, introducing a rate of relaxation towards the local mean direction, under 
the form of a new parameter z/, which can be viewed as a frequency of interaction 
between a particle and its neighbors. It is therefore possible to derive a kinetic 
mean-field limit of this model. Then they develop a method, defining the notion 
of generalized coUisional invariants, which allows to derive the formal limit of this 
kinetic mean-field model, at large scale in space and time. This continuum limit 
is a non-conservative system of PDE for the local mass and the local orientation. 
Moreover, this system is proved to be hyperbolic. 

The goal of the present paper is to confirm the ability of this type of macroscopic 
model to describe the large scale dynamics of systems of self-propelled particles 
with orientation alignment, and to show that the notion of generalized collisional 
invariants is well adapted to derive this model from the microscopic mechanism 
of alignment. This was shown in [17] for a different type of alignment, based on 
the curvature control (for a model of displacement introduced in [16], designed to 
fit biological experiments [23]): the method which uses the generalized collisional 
invariants is successful to derive a macroscopic model which is the same as the 
"Vicsek hydrodynamics" of [15], but for the definition of the coefficients in the model. 
We will show that this is also the case when we slightly modify the individual model, 
in order to be more coherent with some numerical observations, and to model the 
infiuence of an angle of vision. 
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One of the properties of the macroscopic model of [15] which fails to represent the 
numerical observations is that the local equilibria have a constant order parameter. 
This order parameter is indeed only related to the ratio between the frequency p of 
interaction and the intensity d of the noise. In numerical experiments on the Vicsek 
model [11], the order parameter depends on the local density of particles: one can 
observe, at large time, formation of travelling bands of high density, strongly ordered, 
moving through a disordered area of low density. 

The first refinement on the model will be to define the local density p in the 
particular model as the mean number of particles in a neighboring area, and to 
make the parameters v and d depend on p. In a modelling point of view, this could 
be interpreted as the fact that, due to social pressure for example, it is more likely to 
move and update the direction when there is a large number of particles around (so p 
increases with p), and that the fluctuation in the estimation of the mean velocity 
is smaller when a lot of particles are taken in account in the neighborhood (so d 
decreases with p, see for example the way the vectorial noise is defined in [24, 11]). 
This dependence on local parameters for the noise parameter has been introduced 
in other models of collective behavior [19, 41]. 

In [15], the parameter v does also depend on the angle between one particle 
direction and its target direction. For convenience, we will not take this in account 
for most results but, in some simple cases, computations have been done and will 
be given in appendix. 

Another property of the mean-field model of [15] is that the type of local equilib- 
rium for the rescaled model is unique. The loss of this uniqueness could play a role 
to understand the formation of patterns such as the travelling bands [20], and in 
some models of rod-like particles, we have indeed bistability in a certain regime [31]. 
This is not the case here, and we will see in Section 3 that we still have a unique 
kind of equilibrium associated to a local density p and a local orientation f2. 

The second refinement is to take into account the influence of an "angle of vision" 
in the model. In the original Vicsek model, the target orientation for a given particle 
is chosen to be the mean orientation of the neighbors located in a ball centered on 
this particle. We will use here a more general kernel of observation which can be non- 
isotropic. This refinement has been proposed in various models of swarming [1, 30]. 

The main result of this paper is that the formal macroscopic limit of this model 
take the same form as the previous "Vicsek Hydrodynamic model" of [15], consisting 
in a conservation equation for the local density of particles, and an evolution equation 
for the mean orientation, which is not conservative (the velocity is constrained to be 
on the unit sphere). This system of PDE is quite similar to the Euler equations for 
gas dynamics, but presents some specific issues, for example there are two different 
velocities of propagation. The difference between our model and the macroscopic 
model of [15] relies on the definition of the coefficients of this model, and on the fact 
that they depend on the local density. This last feature allows the model to lose the 
property of hyperbolicity. 

In Section 2, we present the individual model and the final macroscopic model, 
focusing on how the two refinements are taken in account, and what are the conse- 
quences at the macroscopic level. 

In Section 3, we provide elements of the derivation of this macroscopic model. 
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following the method of [15], but emphasizing the details which are specific to this 
study. We also give the method in a general n-dimensional framework (previously 
the method was only done in three dimensions). The case of the dimension 2 is 
special, since we are able to give an explicit expression of the coefficients. 

In Section 4, we study the properties of the macroscopic model. We prove that 
when one of the coefficients is negative, then the system is not hyperbolic. We 
describe the region of hyperbolicity for a system which depends only on one space 
variable, and we discuss the influence of "angle of vision". 

Finally, in Section 5, we provide a general method which gives an asymptotic 
expansion of the coefficients in any dimension, in the limit of a small or a large 
concentration parameter. With these expansions, we are able to study the qualitative 
influence of the "angle of vision" in the flnal macroscopic model, and to give examples 
for which the hyperbolicity is indeed lost. 

2 Individual and continuum dynamics of a modi- 
fied Vicsek model 

We start by presenting the individual-based model and the continuum model we 
obtain in the limit of a large number of particles, when observed at large scale, in 
space and time. Elements of the derivation of this macroscopic model will be given 
in section 3. 

2.1 Starting point: particle dynamics 

Here, we briefly recall the time-continuous version of the Vicsek model, and introduce 
how we take into account the anisotropy of observation and the dependence on the 
local density for the rate of relaxation and the intensity of the noise. 

We consider a system of particles with positions Xk in (with k E |1, A^]) 
and orientations in the unit sphere Sn-i; which we will simply write §. 

For each particle, we flrst deflne a local mean orientation (considered as a 
target direction) and a local density pk- In the original model of Vicsek [40], the 
mean orientation ujk is computed on all the neighbors within a given radius R. Here 
we take the mean according to a kernel of observation K, which can be more general 
than the indicator function of the ball of radius i?, as in the time-continuous version 
of [15]. The kernel therein depends only on the distance between the given particle 
and a given neighbor, so the reflnement here is that it can also depend on (the cosine 
of) the angle between the orientation of the flrst particle and the right line joining 
the two particles: 

J - 1 ^ 

= TT^, where Jk = i\^J " iv'Ivj ' ^k) (^j- (2.1) 

\Jk\ .^^ ^ \ J k\ y 

For example, to take into account only the neighbors located "in front", and 
within a given radius R, of one particle, the kernel would be K{r, 7) = l^riiR}^{-y^o}- 
We proceed in an analogous way to compute the local density pk, which may use 
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another kernel K: 

1 ^ ~ 

~P^ = J^Y.K - X,l ■ uu) . (2.2) 

We now turn to the dynamics of the particle system. The k^^ particle moves at 
constant speed 1, following its orientation Uk- This last one relax towards the mean 
orientation of its neighbors, with rate p (depending in the local mean density pk), 
under the constraint that Uk is of norm 1. Finally, this orientation Uk is subjected 
to a Brownian motion (see [26] for more details on how to define such an object on 
a Riemannian manifold, such as the unit sphere here) of intensity d, which will also 
depend on the density pk- The model takes then the form of 2N coupled stochastic 
differential equations, which have to be understood in the Stratonovich sense: 

dXk = Ukdt, (2.3) 

duk = u{pk){ld - Uk ® u]k)uJk dt + \j2d{pk) (Id - Uk ® Uk) o dB^, (2.4) 

where (B^) are independent standard Brownian motions on M". 

Here we denote by Id — ® ujk the projection on the plane orthogonal to Uk, 
that is to say (Id — a; ® uj)v = v — {v ■ u)u. This projection is necessary to keep Uk 
on the unit sphere. The term (Id — a;^ ® Uk)oJk can also be written Va;(w ■ C0k)\uj=uikj 
where V^; is the tangential gradient on the unit sphere. The deterministic part of 
the SDE (2.4) can then be written — i'(^pk)W ui{uj • uJk)\uj=LLJki which is indeed a 
relaxation towards Uk (where the function uj ^ u ■ Uk reaches its maximum), with 
rate u{pk). 

Since the local density p only appears through the functions ly and d, we can 
assume the following normalization for the kernel K: 

[ Kin^-uj)d^ = l. (2.5) 

This normalization condition (which does not depend on G S) means that the 
density is chosen to be 1 in the limit of a uniform distribution of the N particles in 
a region of unit volume. This is not necessary to take a similar condition for the 
kernel K, since Uk, defined at equation (2.1), is independent of such a normalization. 

In [15], the relaxation coefficient u depends on (the cosine of) the angle between 
the orientation of one particle and the target direction, in order to take into ac- 
count the "ability to turn". This would amount to replace iy{pk) by v{pk,0Jk ■ ^^k) 
in (2.4). With our new features here, this would involve many more computations, 
but, following exactly the same method, this leads to the same conclusion. For sim- 
plicity here, we will work without this dependence. We will only present the final 
results with this dependence in some special cases, and add remarks to explain the 
difference in the steps of the derivation of the macroscopic model. 

Numerical simulations tend to show that this time-continuous individual based 
model present the same behavior at large scale as the discrete one (for example the 
formation of bands, as in [11]), in the case where z/ and d are constant, and the ob- 
servation kernel is isotropic, as in [15]. We can expect to observe the same behavior, 
even when u and d depend on the local density p. More precise investigations on the 
numerical comparison between the original discrete and the present time-continuous 
dynamical systems are in progress. 
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2.2 The continuum model 



In this paper, following the approach of [15], we derive, from the particle dynam- 
ics (2.3)-(2.4) introduced in the previous subsection, the following continuum model, 
the functions p{x,t) > and fl{x,t) e § describing the average density and particle 
direction at a given point x G M": 



where the functions Ci, C2, and A will be specified later on: see (3.17) and (3.19)- 



This system of first order partial differential equations shows similarities with the 
Euler system of isothermal compressible gases, but also some important differences. 

Equation (2.6) is the conservation of mass: the density p moves through direc- 
tion Q with velocity Ci(p). We will see that this velocity, taking values between 
and 1, plays the role of an order parameter: when the directions of the particles 
are strongly correlated (close to fl), density moves with velocity close to 1. This 
order parameter depends only on the ratio between i^(p) (the alignment strength) 
and d{p) (the noise intensity). 

Equation (2.7) describes the evolution of the direction Q, the norm of which is 
constrained to be constant (the projection Id — f2cS)f2 insures that the dynamics take 
place on the hyperplane orthogonal to Q). This constraint implies that the equation 
is not conservative which is the counterpart of the fact that, at the microscopic 
level, the only conservative quantity is the mass. The perturbations of this vector 
travel with velocity C2(p), influenced by a term playing the role of pressure due to 
the density, of intensity A(p). It is important to see that in general (and contrary to 
the classical Euler system), the two convection speeds Ci and C2 are different, which 
means that the perturbations on the mean orientation do not travel at the same 
velocity as the "fluid". 

This macroscopic model is the same as the "Vicsek hydrodynamics" of [15], 
except for the definitions of these speeds Ci and C2, and of the parameter A. This 
confirms the ability of this model to describe the global dynamics of systems of self- 
propelled particles with constant speed and alignment interactions (this was also the 
result of [17]). 

However, the parameters depend here on the density p, and their expressions 
are slightly different (due to this dependence and to the anisotropy of the kernel of 
observation). This leads to strong differences in their behavior, as we will see in 
Section 4, devoted to the investigation of the properties of (2.6)-(2.7). For example 
the parameter A can be negative, which implies the loss of hyperbolicity. And 
because of the non isotropy of the observation kernel, the convection speed C2 can 
take a large range of values, from negative if the kernel is strongly directed forward, 
to higher than Ci, if particles are more influenced by neighbors behind them than 
those in front of them (this has been observed in locust migratory bands [5] , where 
the individuals have "cannibalistic interactions" and avoid to be eaten by those 
approaching from behind). 



dtp + (Clip) pQ) =0, 

p {dtn + C2(p)(fi ■ + A(p) (Id - ® n)v^p = o 



(2.6) 
(2.7) 



(3.20). 



6 



3 Elements of the derivation of the continuum 
model 



This derivation proceeds like in [15] but there are significant differences due to the 
additional complexity. In this section, we briefly recall the method of Degond and 
Motsch and will focus on the points which are specific to the present study. 

The derivation proceeds in several steps. The first one consists in writing a 
kinetic version of the particle dynamics. 

3.1 Step 1: mean-field model 

Let f{x,u,t) be the probability density of finding one particle at position x G W^, 
orientation oj E E> and time t ^ 0. The mean-field version of (2.3)-(2.4) is given by 

dtf + u- VJ + V. ■ (F/) = V. ■ (v/^V^V^/)' (3-1) 

with 

F{x, u, t) = v{p) (Id — w ® u)uj{x,u, t), 

x\, ]^ ■ f{y,v,t)dydv, 



2;|i|zf| -w) V f{y,v,t)dydv . 

The first equation (3.1) is the so called Kolmogorov-Fokker-Planck equation. The 
force term F{x,u,t) corresponds to the orientation interaction. 

First of all, if there is no noise (that is d{p) = in equations (2.3)-(2.4), which 
become ordinary differential equations), the formal derivation of this system is easy: 
the usual methodology shows that the empirical distribution (see [38]) satisfies the 
equation (3.1), with d = 0. 

Difficulties appear when the noise is added. A method consisting in writing the 
BBGKY hierarchy (see [25], applied to the Cucker-Smale model of self-propelled 
particles) would not in that case reduce to a evolution equation involving only the 
one-particle and the two-particles distributions, since the interaction is not a sum 
of binary interactions. 

We could slightly change our model to make the interaction as a sum of binary 
interactions, replacing uik by Jk (defined in equation (2.1)) in the system of parti- 
cles (2.3)-(2.4). In that case the model present an interesting phenomenon of phase 
transition and is the subject of current work in collaboration with P. Degond and 
J.-G. Liu [14] (the homogeneous version has been studied with him in [21]). In that 
case, writing the BBGKY hierarchy and using exchangeability of particles gives 
a system of evolution equations, the first one involving only the one-particle and 
the two-particles distributions. The classical assumption of propagation of chaos 
amounts to consider the two-particles density as the tensor of the one-particle den- 
sity / by itself (in the limit of a large number of particles, two particles behave as 



p{x,u,t) 
u(x, u, t) 
J{x, u, t) 



K{\y 



J{x, u, t) 



\J{x,UJ,t)\' 



K{\y- 
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if they were independent), and this gives exactly the evolution equation (3.1). Ac- 
tually, it has been recently proved [9] that the mean-field limit of this model is the 
partial differential equation (3.1), where u is replaced by J in the definition of the 
force F (in the case where u and d are constant). The main point to derive this limit 
is to adapt the classical theory of propagation of chaos [32, 39] in a framework of 
stochastic analysis in a Riemannian manifold (the unit sphere in the present case). 

In the case of a non-linear diffusion, some results are given in [10] for other sys- 
tems of self-propelled particles, under assumptions which would have to be adapted 
in our framework. We can expect to have conditions such as to be Lipschitz for 
the function ^/d, and to be Lipschitz and bounded for the kernel K. Since we 
use the Stratonovich formulation in order to work on the unit sphere, we get the 
term Vo; ■ {\J d{p)V u}\J d{p) f) (instead of Ay{d{p)f) when the velocity w e M"^ satisfies 
the SDE in the usual Ito formulation, see the sections 4.3.5 and 4.3.6 of [22] for the 
correspondence) . 

Finally, when the drift is not under the average form, it is sometimes possible 
to get a mean-field limit, under regularity assumptions on the coefficients [35], or 
with weaker assumptions, but assuming uniqueness of the solution of the mean- 
field model [34]. These results could to be adapted in the framework of stochastic 
differential equations on the sphere, but dealing with the singularity of u (when J 
is close to zero) seems to be slightly more complicated. 

With these considerations in mind, it is however very reasonable that the limit 
of the particle system (2.3)-(2.4), when the number of particles is large, is given by 
the mean-field model (3.1). So we start with this model as a base for the derivation 
of the continuum model. A rigorous proof of the derivation of such a mean-field 
model from the individual dynamics is left to future work. 

Remark 3.1. If we want to take into account some "ability to turn", we just have 
to replace u[p) by iy{p, u ■ u). 

The next step consists in observing this system at large scale, in both space in 
time. 

3.2 Step 2: hydrodynamic scaling 

The hydrodynamic scaling consists in the same rescaling for the time and space 
variable. We introduce a small parameter e and we set x' = ex, and t' = et. 
We define /^(x',a;,t') = f{x,u,t), and we rewrite the equation (3.1) in this new 
coordinates. 

The kinetic equation has the same form, with a factor e in front of each of the 
terms with space or time derivative: 

eid.r + CO ■ v^^n + ■ {F^n = ■ {,id{f)v^,ld{f)n, 

with 

F'{x\uj,t') = v{p^) {Id- u®u)uj'{x\uj,t'), 



8 



where the local rescaled density and orientation are given by 

nx\ to = l^^^^^^^ K ^ . . r d. . 

The important point is to realize that the average density fF and orientation Uj^ 
now depend on e, and can be easily expanded in terms of e, the non-locality only 
appearing at high order. Omitting the primes for simplicity, we have the following 
expansions, the proofs of which are given in Lemma 3 of Appendix A.l: 

cj, t) = n^ix, t) + ea{u- V^) n'{x, t) + 0{e^) , 
p^(a;, u, t) = p^(a:, t) + eau ■ V^p^(x, t) + O(e^) , 

where = p/s and Q"^ = Qfe are the local density and mean orientation associated 
to the function (these quantities, depending only on the space and time variables, 
are related to the first moments with respect to the variable u) given by 



pj= f{.,io)duj, (3.2) 



^/ = "TH"! with J/ = / uf{.,uj)doj, (3.3) 
\jf\ Jues 

and the constants a and a depend only on the observation kernels K and K. These 
constants are positive if the kernel is directed forward, and the more acute the "angle 
of vision", the bigger the constant related to the kernel. 

Now we can introduce these expansions in the mean-field model, and after some 
easy algebra, the rescaled model can be written in the form 

eidtf + CO ■ V.r + a Pin + 5 Pif)) = Qif) + Oie') , (3.4) 

where Q, P and P are the operators given by the following equations (where z> and d 
are the derivatives of u and d with respect to p): 

Qif) = - ^(p/) V. ■ ((Id -u^ u)nff) + d{pf)A^f, 

P{f) = u{pf)V^ ■ ((Id ~u®u){{u- V,) fi/)/), 
P{f) = u{pf)V^ ■ {{u ■ V^pf) {\d-u® uj)Qff) 

- dipf)\/^ ■ (i(Id -u® u)\/^pff + (u ■ W^pf) V^/). 

Notice that the operator Q (giving the only term of order in e) only acts on the 
variable u, and the study of its properties will be important for the following. 

Remark 3.2. Ifu also depends on u-u, the expression of the operator Q is the same 
with v{pf,uj • Vtf) instead of v{pf). But then the expressions of P and P complicate 
in a significant way, since there are also terms with the derivative of v with respect 
to this second variable. 

Now we are ready to study this system when £ — > 0. 
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3.3 Step 3: limit as e — 

This is the main step, where we give the hnk between the continuum hmit (2.6)-(2.7) 
and the rescaled kinetic equation (3.4) of the particle dynamics. 

Theorem 1. The limit when 5 — )■ o/ is given (formally) by p = pM^(^p-^^ 
where p = p{x,t) > is the total mass of p and Vt = Q{x,t) G § its mean orienta- 
tion: 

p{x,t) = / f°{x,u,t) dw, 

^=-n", j{x,t) = f'^{x,uj,t)ujduj, 
131 -'^es 

where M^^ is a given function of u ■ Q and ^ = ^ which will be specified later 
on (see (3.5)^. Furthermore, p{x,t) andQ{x,t) satisfy the following system of first 
order partial differential equations: 

dtp + V, ■ {cipQ) = 0. 
p {dtfl + C2(fi ■ + A (Id - 1] ® ft)V^p = 0, 

where the convection speeds Ci, C2 and the parameter A depend on p. Their expres- 
sions will be given in this section (see (3.17) and (3.19)-(3.20)j. 

The method to obtain this result follows closely [15], and is only summarized 
here. We will focus on the details which are specific to this study. 



3.3.1 Equilibria 

The first important point is to determine the null space S of Q, since it contains the 
limits of (3.4). We find, as in [15], that it is a n-dimensional manifold consisting of 
functions analogous to Maxwellian distributions in the classical Boltzmann theory: 

S = {pM^^p)n{uj)\p>0,neS}, 

where 

= ^ > and M.„(.) = (3.5) 

The main difference with [15] is the dependence on p for this equilibrium in a nonlin- 
ear way, coming from the dependence of v and d on p. This will result in additional 
terms in the computations, and so in additional terms in the expressions of the 
constants in the macroscopic model. 

The normalization constant f^e*^^'^ du depends only on k (not on Vt) and so 
the total mass of Mi^q{uj) is 1 and its mean direction is 1), that is to say Pa/^^ ~ 1 
and i^M^n — ^- Indeed we can easily compute the flux ju^n ^^^i^ equilibrium, 
defined by (3.3), and we get: 

jM^o = (cos^)a/«^^, (3.6) 
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where for any function 7(cos6'), the notation {'j{cos6))m^ stands for the mean of 
the function oj 7(0; • fl) against the density M^n, i-e. 

(7(cos^))m. = / M.oM7(a;-fi)du;- ^^^^'"■^^'''"■''^'^ 



je§ /g e'^'^'*^ do; 

Notice that {j{cos9)) depends only on k, not on Q: 

/o-7(cosg)e--^sin"-^gd^ 
(7(cos^))m. = j^^e-co..,,,n-.gd^ ■ (3.7) 

Remark 3.3. In the case where v depends on p and u -Q, we have to replace in all 
this point kuj ■ by to ■ il), where k{p, p) = Jq ^^^dr. 



3.3.2 Collisional invariants 

The second important point is the determination of generahzed colhsional invariants. 
Indeed, since there is no other conservation relation than the conservation of mass, 
the collision invariants reduce to the constants, and the integration of the equation 
against these invariants only gives one equation, which is not sufficient to describe 
the behavior of the equilibrium (which lives on a n-dimensional manifold). The main 
idea in [15] was to overcome this problem with a generalization of the concept of 
collisional invariants. 

A collision invariant is a function ■0 such that for all function / of u, the inte- 
gration of Q{f) against "0 is zero. So we ask for a generalized invariant to satisfy 
this definition only for a restricted subset of functions /. In the case where the 
dependence on p in the equilibria is linear, restricting to all functions with a given 
orientation Q is sufficient to obtain the remaining equation. Here we also have to 
restrict to functions with a given density too (actually, we impose a given k{p))- We 
will have then a set of generalized coefficients indexed by G S and k > 0. 

More precisely, to have a good definition, we have to work with linear opera- 
tors (this point has been mentioned in [17], since the result given in [15], with the 
definition therein, was slightly incorrect). We first define the linear operator L^n by 

L^nif) = -A^/ + kV^ ■ ((Id - a; ® io)nf) = -V. ■ 

and then the generalized collisional invariants C^n (associated to k G M and G §) 
as the following vector space: 

C^n = (^1 / L^aif) ijdu = 0, V/ such that (Id - ® n)jf = 

We remark that the operator Q{f) can be written as Q{f) = —d{pf)Lf^(^pj)Q_j{f). 
Hence, for any generalized collisional invariant if) G C^f^, we have 

V/ such that Vtf = Vl and K,{pf) = k,, Qif) i^du = 0, (3.8) 

and this is the only property of generalized collisional invariants we will need in the 
following. 

The computation of the set of generalized collisional invariants has been done 
in [15] in dimension 3. We give here the general result in any dimension. 
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Proposition 3.1. Structure of the generalized collisional invariants. 

Any generalized collisional invariant ip associated to k E M. and fl E S has the 
following form: 

tjj{u) = C + h^{u ■ Vl)A ■ oj, 

where C E is a constant, the vector A G MJ^ is orthogonal to Q, and /i^ is a given 
positive function on (—1,1), depending on the parameter k, which will he specified 
later on. In particular, the generalized collisional invariants form a vector space of 
dimension n. 

Proof. We first reAvrite the set {/| (Id — 17 (g) fi) = 0} as the set of functions / such 
that for all A G M" with A ■ f2 = 0, we have that j^A ■ ufdoj = 0. Finally this is 
the orthogonal of the set {u ^ A ■ u for A ■ Q = 0}, for the usual inner product 
on L^(S). We can then rewrite the set of generalized collisional invariants: 

C.n = l^P\ I f Ll^^duj = f e {uj ^ A - uj ioi A - ^ = Q}"-] 

= ^e{{uj^A-ujioi A-n = 0}^)^} 

= {^iL^f^ ^(w) = A-uj with A ■ = 0}, 

the operator being the adjoint of the operator L^q, which can be written 

Ll^ V = -^.>^ - ■ Vc.^A = -T^V^ • (M^nV^tP). (3.9) 

It is then easy to show that the problem L*^ '^/'(w) = A ■ u, for A ■ Q = has a 
unique solution in the space if^(S) (functions of H^(E>) with mean zero), using Lax- 
Milgram theorem and the Poincare inequality. Hence, if we show that this solution 
has the form iplu) = hi^{Q ■ u))A ■ u, the solutions in the space H^{Ei) are equal to 
this solution plus a constant C. 

We search a solution of this form. We identify Q with the last element of an 
orthogonal basis of W^, and S„,_2 with the elements on the unit sphere § which 
are orthogonal to Q. We can then write u = cos6Q + sin6'f, where v G Sn-2 
and 6 G [0, tt] (this decomposition is unique when u is different from Q and —Q). In 
this framework we try to find a solution of the form ip{uj) = /ik(cos^) sin^y4 ■ v. 

For ip{uj) = g{6)Z{v), we have, in dimension n ^ 3: 

V.ipico) = g'{9)eeZ{v) + ^V,Z{v), 

sm o 

where the unit vector eg is Vo;^ = — ^j^(Id — w ® u)Q. A tangent vector field can 
always be written acg + A where ^ is a vector field tangent to the sphere Sn-2, and 
we have 

Vo, ■{aee + A) = sin^"" 9 ^^(sin"-^ 60) + ■ A . 

sm U 

Finally we get, using the second part of (3.9), 
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In our case we have Z{v) = A ■ v, so we get A^Z = —{n — 2)Z (this is a spherical 
harmonic of degree 1 on §n-2)- So L^.^'ip is also of the form L*^g{9)Z{v), where 

LlgiO) = -sin'-^ ee-'^^"^'^,{sm^-' 9e^^°^'g'{e)) + ^^(^). (3.10) 

Finally, solving {h{u ■ Q)A ■ u) = A ■ u comes down to solving 

L;^ = sin^, with g{9) = h{cos9) sinO. (3.11) 

When y4 7^ 0, it easy to see that the function uj i— t- h{u-Q)A-u belongs to H^{S) if and 
only if the function g : 6 ^ h{cos6) sin 6' belongs to the space V (a "weighted -f^o") 
defined by 

V = {g\{n- 2)(sin^)t-2^ e L\0, vr), {sm9)^-^g e H^,{0, tt)}. (3.12) 

Using again Lax-Milgram theorem in this space V, we get that the problem (3.11) 
has a unique solution, denoted g^^, which is positive (by the maximum principle). 
Writing h^in) = (1 — /i^)~2 (7K(arccos(//)) gives that = hi^{u ■ Q) A ■ u is a, 

solution to the partial differential equation L*^ V^(^) = A ■ u. We could write the 
elliptic equation on (—1,1) satisfied by /i^ to have another definition, but this does 
not give a more elegant formulation. 

In the case of dimension 2, we write iIj{u) = g{9)A ■ vq, where Vq is one of the 
two unit vectors orthogonal to fl and g is an odd 27r-periodic function in Hl^^{R), 
which can be identified with a function g G i^fo(0,7r) = V. We still have that the 
elliptic problem L*^^tjj{uj) = A ■ u is equivalent to (3.11) with g E V, with the 
same definitions (3. 10)- (3. 12) of and V . But since this elliptic equation reduces 
to {e'^'^"^^ g'{0)y = — sin 6 e'^ ^ , we now have the following explicit expression of g^: 

9 vrj^-!!!!^ r3i3i 

□ 

Remark 3.4. If we take into account the "ability to turn", we just replace kcos9 
in equation (3.10) by k{cos9). In dimension 2, we still have an explicit expression: 



where 



g^{9) = - / e^(^°^'^)-^(™'^^)sin0d(/)dv?, (3.15) 

Jo J ip 

gf(9) = Te-^^^'^^Mv?. (3.16) 



3.3.3 Computation of the limit as £ — )• 

The third and final important point is taking the limit e — )■ in the equation (3.4), 
after integrating against the collision invariants. Since we do not have results of 
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existence, uniqueness and regularity of the solution, all the limits in this section 
have to be understood as formal limits. A rigorous proof of convergence is left to 
future work. 

When £ —J- 0, if we fix x and t > 0, we have that Qif^), as a function of u, 
tends formally to zero, so tends to an equilibrium of the operator Q, of the 
form pM^(^p-^^, where p > and f2 G S are given functions of x and t. This is the 
first part of Theorem 1. So we have — )• p, and Q'^ ^ Q. When there is no possible 
confusion, we will write k for K,{p). 

For the mass equation, we use the constant invariant: we have, since the opera- 
tors Q, P and P are given as the divergence (with respect to w) of a function, 

/ Qif) duj= [ Pif) du= [ Pif) du = 0. 

Juj& Juj& Juj& 

Hence, integrating the equation (3.4) with respect to u, we get: 

dtp' + V.-f = 0{e). 

Actually, we can even replace the 0{e) by zero in this equation since in the original 
model (3.1) we have conservation of mass. We get in the £ — )■ limit: 

dtp + \/^-iciin)pn) = 0, 

where (see (3.6)): 

Jo 

This gives the second part of Theorem 1, with the equation on p and the definition 
of Ci. 

To get the equation on Q, we use the non-constant part of the collisional invari- 
ants. By Proposition 3.1 and using the result at equation (3.8), we get that for all A 
such that A ■ = 0, we have 

/ Qir)K(p^){u-n'')A-ujduj = 0. 

So we have that the vector = Xjg§ Q(/^)^K(pe)(^ ■ fi^)a;dw is orthogonal to A 
for all A orthogonal to Q, that is to say that X'^ is in the direction of fi^, which is 
equivalent to (Id — fi^) = 0. Using (3.4), we get that 

X' = e [ {dtf + uj ■ V.r + aP{f) + 5P(r)) K(p^){u ■ fi^) ujdco + 0{e^). 

Dividing by e and taking the limit e — 0, we get (Id — Q ® Q) X = 0, where 

X = f {dt{pM,n)+uj-V,ipM,n)+aP{pM,n)+aP{pM^n)) K{uJ-n)ujdco . (3.18) 

The main point is then to compute (Id — Q ^ Q) X, in terms of p, Q and their 
derivatives, using mainly the chain rule. The computation is similar to [15] for some 
terms, but additional work is required for the terms coming from the nonlinearity 
of Mf^n in p and the operators P and P. We give the result of the computations 
under the form of a proposition: 
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Proposition 3.2. (Id — ® f^) = 0, where X is given in (3.18), is equivalent to 

p {dtn + C2{n ■ v^)i]) + A (Id - ® n)v,^p = o, 

where 

C2 = Cl — a d {nci + K. {cos^ 9) ) , with ci = {cos9)^^ , (3.19) 
A = i + [ci - Cl + «(/(/€ (sin^^)^^ - nci)] + pd{n - 1 + KCi) , (3.20) 

with the notation 

/o^7(cos^)/i«(cos^)e^™^^ sin"^d^ 



(7(cos^)) 



This proposition is exactly the last part of Theorem 1, with a precise definition 
for coefficients C2 and A, and this ends the derivation of the continuum model (2.6)- 
(2.7). 

The computations to get this result are given in Appendix A. 2, the idea is to 
compute (Id — Q ^Q) X using the chain rule and the change of variable u ^ {6, v) 
where u = cosOQ + sin^v , with v orthogonal to Q, which simplifies a lot of terms. 

Remark 3.5. The computations have also been done in the case where u depends 
on u ■ Q (and not on p) and where d is a constant. We get the same results, except 
that the constants are given (with analogous definitions) by: 

ci=(cos^)^,^ (3.21) 

C2 =(cos6') — a(z/cos^ 6 — u' cos 6 sin^ 6) 

- ad (n COS 9 + '^{{n + 2) cos^ 9-1) - ^cos^sin^^J 



(3.22) 



A=d(i)^^, (3.23) 

where here we use the notation 

£ 7(cos 9)u{cos 9)h-{cos 9)6^^"''^^^ sin" 9 d9 



(7(cos 9)) 



A/- 



/J'//(cos6')/i;j(cos6')e'^(=°'^^) sin^^d^ 



Since v is supposed to be positive, the constant A is positive, and we will see in the 
next section that its possible change of sign with the dependence on p is important. 
This is why we focus on the dependence on p and not in u ■ fl in this article. 

4 Properties of the macroscopic model 
4.1 Hyperbolicity 

We recall here the macroscopic model (2.6)-(2.7): 

9tP + V,-(ci(p)pfi) =0, 
p {dtil + C2(p)(fi ■ V^)fi) + A(p) (Id - (g) fi)V^p = 0, 
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where the functions Ci, C2, and A are given by (3.17) and (3.19)-(3.20). A first remark 
is that it is not possible to do another scahng to get rid of ci, hke in [15], because ci 
depends on p. 

The main result about this model is that li d 01 v depends on p, the coefficient A 
can become negative in some regions of the state space, and in that case the system 
loses hyperbolicity. Let us first discuss here the interest and the problems due to 
the non-hyp erbolicity. 

The first thing to remark is that the model is not always well-posed. Indeed, 
in general, we cannot ensure that a solution will stay in the region of hj^erbolicity 
for all time, even with smooth initial conditions in the hyperbolic region (actually, 
even with hyperbolicity everywhere, dealing with the discontinuities is a challenging 
issue, see [33]). 

The property of hyperbolicity is linked with the fact that perturbations propagate 
with finite speed. Here the presence of a region of non-hyperbolicity means that we 
could have propagation with infinite speed across this region. This leads to a second 
remark: it may be possible to construct non-classical shocks, using the crossing of a 
zone of non-hyperbolicity, see [29], and [28]. The interest is that we may construct 
some travelling waves, as observed in [11]. Actually we did not manage to construct 
such solutions yet, this is part of our future work. 

We should also construct models with formation of coherent structures from 
such non-hyperbolic models, if we could use stabilization with diffusion. But here 
the expansion at higher order in e in the rescaled mean-field model (3.4), in order 
to obtain diffusion terms in the macroscopic model becomes too much complicated 
to perform some study (see [18] for the case of the original model of [15]). 

We now turn to the description of the regions of non-hyperbolicity. We consider 
a system satisfying (2.6)-(2.7), but evolving only along one space direction G S 
(the density p and the orientation f2 depending only on 2; = ■ x and t). We 
write then Vt = cosO + sin ^-y, where v G S„_2 (identified to the set of unit vectors 
orthogonal to e^). In this framework, the system is equivalent to 



In the special case of dimension 2, the system reduces to (4.1)-(4.2), with 6 G (— vr, tt) 
and VL = cosd Cz + sin^vo, where vq is one of the two unit vectors orthogonal to Cz- 
The general definition of a quasilinear hyperbolic system [37] gives that the 
system (2.6)-(2.7) is hyperbolic if and only if this system (4.1)-(4.3) is hyperbolic 
for all unit vector G S. We give the result in the following statement: 

Theorem 2. Hyperbolicity. 

• The system (2.6)-(2.7) is hyperbolic if and only if \{p) > 0. 

• The system (4.1) -(4.3) is hyperbolic if and only if 



dtP+ dz{pci{p) cos 9) = 0. 
p[dticose) + C2{p) cos (cos ^)] + A (1 - cos^ 9) dzp = 0. 
dtv + C2(p) cos9dzV = 0, with |t>| = 1 and Cz ■ v = 0. 



(4.1) 
(4.2) 
(4.3) 



A(p) > or 




ifX<0, 
ifX = 0. 



(4.4) 
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where Oi{p) = ^(pci(p)) = ci(p) + pk (^{cos'^9)m^ - (cos6')^^J. 



Proof. The system (4.1)-(4.3) can be written as the following first order quasilinear 
system of partial differential equations 





with 



+ A{p, cos 9, V 



( C3{p) cos 9 ci{p)p 
-sin^6' C2(p)cos6' 



0, 



v4(p, cos6', v) 



\ 





V 



















C2(p) cos6'Id„_2 



and this system is hyperbolic in case A > 0. The eigenvalues are 7± and 70 (of 
multiplicity n — 2), given by 



7o = C2 COS 9, 7± = 2 



1/2' 

(c2 + C3) COS 6* ± { (c2 - cg)^ cos^ 9 + 4Aci sin^ 9^ 



Now if A < 0, asking 7-1- to be real and distinct is exactly equivalent to the equa- 
tion (4.4). In this case the matrix A is diagonalizable. If 7_|_ = 7„, then A is 
diagonalizable only if its top left corner 2x2 submatrix is scalar (only one eigen- 
value), which is not the case since ci{p)p > 0. For the same reason, if A = 0, we 
immediately get that A is diagonalizable if and only if the first two diagonal coef- 
ficients C2(p) cos^ and C3(p) cos 6* are different, which ends the proof of the second 
statement. 

Now we turn to the general case. If A > 0, the system (4.1)-(4.3) is hyperbolic for 
all unit vector G S, which gives that the system (2.6)-(2.7) is hyperbolic. Suppose 
now that the system (2.6)-(2.7) is hyperbolic in an open region of the state space, 
with A ^ at some point (p, Q). Since n ^ 2 we can find such that e^-fi = 0. Then 
we have cos ^ = 0, which gives, by the condition (4.4) that (p, Q) is in the region of 
non-hyp erbolicity of the problem (4.1)- (4. 3), and this is a contradiction. □ 

Actually, the positive functions d and ly being arbitrary, it is possible to have a 
lot of qualitatively different shapes for the region of non-hyperbolicity of the reduced 
system (4.1)- (4. 3). We give here some examples in the case of dimension 2, where the 
coefficients are easy to compute numerically (using the explicit formulation (3.13) 
for g^). 
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h' = p, d = a = 0, a = —0.03 



TV 




Pi 2 A P-2 6 8 



Density p 



Figure 1: The region of hyperbolicity can be of the form (0, pi) U (p2, +00) 




4.2 Influence of the anisotropy 

On the final macroscopic model, the influence of the anisotropy in the observation 
kernels is only visible through the values of the speed C2, and the coefficient A. 

We remark that the parameter a, which is related to the kernel K used to define 
the local orientation a), only appears in the expression (3.19) of the velocity C2, 
making it smaller when a is a large positive constant. The difference between ci 
and C2, which is one of the differences between the macroscopic model (2.6)-(2.7) and 
the classical Euler system, is then enhanced when a is a large positive constant. This 
can be interpreted as follows: if the observation kernel is strongly directed forward, 
then the information on the orientation moves rapidly backward. This could be 
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compared to results on modelling of traffic flows, where the speed of a congested 
phase depends on the distance of anticipation of the drivers (see [4, 8, 13]). 

The parameter a, related to the kernel K which is used to define the local 
density p, only appears in A, and obviously has an influence only if the relaxation 
frequency u or the noise intensity d depends on this density (in the expression (3.20), 
we must have k ^ or d ^ 0). So the anisotropy of the kernel K can have an 
impact on the region of non-hyperbolicity for the system (2.6)-(2.7). The anisotropy 
of the kernel K does not play a role in this global hyperbolicity, but, through the 
condition (4.4), it can change the shape of the region of non-hyperbolicity for the 
one-dimensional reduction (4.1)-(4.3). 

Since the expression (3.20) involves a lot of terms which can take different signs, it 
is not easy to directly quantify the influence of the parameter a on the coefficient A, 
as it was the case for a and C2. In the next section, we perform an asymptotic 
study of the coefficients of the macroscopic model (2.6)-(2.7), as the concentration 
parameter k tends to infinity (in the case of strong alignment, or low noise) or to 
zero (when the noise is high, or the frequency of alignment). 



5 Asymptotic study of the coefficients 

We want to obtain an asymptotic expansion of ci, C2 and A given by the expres- 
sions (3.17), (3.19), (3.20) as the parameter k tends to infinity or to zero. 

Since we do not know explicitly the dependence on p for the coefficients u 
and d, the only quantities we can study in the expressions of this coefficients are 
the averages ci = (cos^)a/«, ci = {cos9)^j , and {cos^ 9) (since the last average 
is (sin^ 9) ^ = 1 — (cos^ ^) m )■ The purpose of this section is to give a method to 
get the Taylor expansion up to any order in k or ^ of the following averages: 

/o^/(^)e'''°'^ sin"-2^d^ 

Jo fi9)Kicos9)e^^'''^ sm''9d9 
/o''/i«(cos^)e'^'=°«^ sin"6'de ' 

where h^, is the function providing the generalized collisional invariants (see Propo- 
sition 3.1). We ffist give the method to obtain the expansion of the ffist type of 
average, and we apply it to get an expansion of Ci in k and -. 



(/(^))m. = 
(/(^))m. = 



5.1 Asymptotics of {f{0)) 



The ffist expansion, when k — 0, is just a basic Taylor expansion. For a function / 
such that /sin"~^^ G L^(0,7r), we define 



h = - r f(9) cos^ 9 sin"-2 9 d9 and a„ = - T cos^ 9 sin"-^ 9 d9. 
pi Jo pi Jo 

Then we get 



s;=obp/^^ + Q(/^^ 
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If we take f{9) = cos 6, we have bp = {p + l)ap+i and integrating by parts, we 
get the following induction relation: {p + 2)(p + n)ap+2 = ttp. Since ai = 0, the odd 
terms vanish and we get 

-1^ + o /io> ^^ + O(k^) 
I 0\ n ' 2n(n+2) \ J i ^ 3 , 5, , - 

Ci = {cOSe)M^ = 1^^^2_^ 0(^4) = n'^ - ^^H^)'^ + 0{f^ )• (5.2) 

We now turn to the expansion of ci when k — oo. We will use the following 
lemma, the proof of which is elementary, see [7] for examples and variants: 

Lemma 1. (Watson's Lemma) 

Let p be a function in L^(0,T), with T > 0, and let Ik{p) = Jo p{t)G~'^^dt. 
Suppose that, in the neighborhood of 0, we have p{t) = t^ {^f^^^ afi + Oit'^^ 
with /3 > — 1. 

Then I^{p) = k"^^^ (e^o^ aiT{l3 + i + + 0{k~^)) as k^oo. 

We use this lemma, after the change of variable t = 1 — cos 6, in the integrals of 
the form [/(^)]« = £ /(^)e'^™'^^ sin'^-^ ^d^. We get 

[fi9)U = r f{arccos{l - t))e^^\2t - t^)^ dt. 

JO 

So if we can expand the function t H- {2t—t'^)^f(arccos{l—t)) in the neighborhood 
of 0, we can apply directly Watson's Lemma to get an expansion of [/(^)]k, and then 
to which gives finally the expansion of (/(6'))m„- 

We take here the example of the function f{6) = 1 —cos 6, so f{arccos{l — t)) = t. 
We want an expansion with two terms (since we have Ci = 1 — {f{d))M^ "we will 
actually get three terms for Ci). We have 

(2t - t^)— = 2—t—{l - It)— = 2—t—{l - ^t + 0{f)) 

Applying directly Watson's Lemma to this function and to the same function mul- 
tiplied by t, we get 

n-3 



[1]. = '^-AS (r(^) - ^r(^)i + o(k-2)) (5.3) 

n — 3 

[/W]. = (r(^) - ^r(^)- + o(«:-2)) . 

K, 2 \ K J 



Since r(p + 1) = pr(p), we finally get 

\fiff\\ y{ ri+\ \ 1 _ n-3 n+1 1 

lf(a\\ — Viy^n'^ _ ^\ 2 I ^ 4 2 K I /nf-,,-3^ 



n — 1 (n — l)(n — 3) 
In particular we get the expansion of Ci as k — oo: 



+ 0(fi:- 



-3^ 



^ ^ _ n-l ^ (n-iK^ ^ (5.4) 



Using this method we can easily get the following lemma, which will be useful in 
the next subsection. 
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Lemma 2. Estimation of {f{6)) m,,- 

Suppose that 9 ^ f{e) sin"-^ 6* belongs to L^(0,7r), and that |/(6')| = 0{e'^^) in 
the neighborhood ofO, with jS > — Then {f{0))M^ = 0{k~^) as k oo. 

We now turn to the method to compute averages of the form {f{9))^^ . 
5.2 Asymptotics of {f{0))^ 

We first decompose /ik(cos6') as a polynomial in k or in whose coefficients are 
polynomials of cos 6 plus a remainder which will be negligible. 

Proposition 5.1. Expansion of h,^. 

We define the two linear operators L and D on the space of polynomials by 

L{P) = -(1 - X^)P" + (n + l)XP' + (n - 1)P 
D{P) = -{1-X^)P' + XP. 

We have the two following expansions: 

N 

/i^(cos^) = Y.Hp{cose)KP + R^q{cos9), 

N 

p=i 

where the Hp (resp. Gp ) are the polynomials of degree p (resp. at most N — p) given 
by the following induction relations (the second one being in the neighborhood ofO): 

Lm = l . /D(Gf)(cos^) = 1 + 0(0 ,,,, 

L{Hp^,) = -D{Hp) \(D(G^^^,) + L(G',^))(cose) = 0(^2(^-^)), ^ ' 

and where the remainders satisfy the following estimations, for any function f such 
that 6 v-f f{6) sina 6 belongs to L'^{0, tt) and such that \ f{0)\ = 0(9'^^) in the neigh- 
borhood ofO: 

{f{e)R^Q{cose)smH)M^ = 0{k^+^) as k^O, 
{f{e)R^,^{cose)smH)M^ = 0{k~^~^-^) as k^oo. 

The proof of this proposition is given in Appendix B.l. The first thing to do 
is to prove that the inductions relations (5.5) make sense to define the sequence of 
polynomials (an induction relation is given in Appendix B.l to compute easily the 
polynomials and Hp). The operators L and D are made so that 

Z*(P(cos^) sin^) = (L(P) + KD{P)){cose) sin 6, 

where the operator L* is defined in (3.10). Since we have L* (/ik(cos6') sin 6*) = sin^ 
by definition, we are then able to obtain the estimates on the remainders, using 
Poincare inequalities in adapted spaces. 
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With this proposition, it is then easy to get an expansion for {f{0))^j , using 
expressions of the form {g{9))Mi^, which can be expanded by the tools of the previous 
section: 



urn 



E%{Hp{cose) sin2^)M,KP 



+ 0{k 



+ 0{ 



as ft; — )■ 0, 



as K — )• oo. 



As an example, we can compute the first polynomials, we get 



Hence, 

(cos^-l)^:2„ 



n — 1 



, Hi 



-X ^2 _ 4 - X 



2n{n- 1) 



Go 



2(n-2) 



;^(cos^sin=^g)M. - 2;^(cos^^sin^g)M. 
f (cos 9{A - cos 6) sin2 e)M^ + (cos 6 sin^ i 



^L((4-cos^) sin2^)M. + ii^(sin 



1 + 



As before, in the second equation, we computed (cos6' — instead of {cos9)j^j 
in order to have a remainder of order 3 instead of 2. 

Finally, we have to compute terms of the form (cos^ 9 sin^ 9)mk- Instead of using 
the method of the previous subsection, we can actually express all these terms in 
function of Ci = {cos9)m^, by integrating by parts. We get 

{sm^9)M. = ^ci, (cos^sin2^)M« = ^(1 - ^ci) 
(cos^^sin^^),,. = (sin2^)M. - (sin^^).,. = ^(ci - ^(1 - ^Ci)). 

Using the previous expansions (5.2) and (5.4), we finally get the expansions of Ci. 



ci = (cos( 



2n(n+2) ^ + ^("^^ 



K 



n+1 _|_ (n+l)(3?i-7) 



2k 



24k^ 



+ 0( 



as K — > 0, 



as K — > oo. 



(5.6) 



In addition, we can compute an expansion of (sin^ 



in order to get expansions for 



the coefficients C2 and A of the macroscopic model, given in equations (3.19)-(3.20). 
Using only Hq and G\ = 1, we get that 



' (cos^ 9 sin^ 



(sm 



(sin^ 9)m. 



+ 0(k) = ^ + 0{k) as k ^ 0, 



(sin^ I 



I (sin^ ^)m. 



0(/€-2) = H+i + 0(/€ 



as K — 7- cxD. 
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So, using the expressions (5.2), (5.4) and (5.6), we get the following expansions: 



Cl 



C2 



2n(n+2) ^ 



-i + 0(k-2)- 



2n+l 
2(n+2) 

azy(l-K-i + 0(/€-2)) 



as K — 0, 
as K — oo, 



I (1 + P ! [-1 + 0(«;^^) + aiy{l + 0(/€^i)) ]) 
+ |ap^z/(l + ^K^i + 0(/€-2)) 

This shows that in any dimension, there are some simple cases were we actually 
have A < and the system loses hyperbolicity, even if 5 = (for example if the 
kernel of observation K is isotropic). For example if we have k = with /3 > 0, we 
have as rho — ^ oo that A(p) = (1 — f3)p~^ + 0(p~^^), which gives that A(p) < if 
we take (3 > 1 and p sufficiently large. 

These expansions also give a more precise estimation on the difference between ci 
and C2 as the noise is small or large: when the kernel of observation is isotropic, we 
have Cl > C2 in the two expansions, in any dimension n. That means that the 
information on the orientation propagates slower than the "fluid". 

Remark 5.1. We can also do an expansion in the more general case where v depends 
on p and u -Vt. When d — > 0, the expansion of the coefficients depends only on the 
local behavior of the function x H- v{p^x) near 1. In Appendix B.2 we give tips to 
perform this expansion. Here we only give the final expansion in the case where v 
and d do not depend on p, so the coefficients are given by (3.21)-(3.23). In this case 
we can suppose u{l) = 1 (up to a rescaling), and denote 7 = y'iX). We finally get, 
when n = 2: 

ci = l- ld + 0{d^), 
C2 = 1 - a + ((1 + f 7)« - l)d + 0{d^), 
A = d+l-fd^ + 0{d^). 



as K — )• 0, 
as K — > 00. 



6 Conclusion 

In this article, we have seen that the introduction of a dependence on the local den- 
sity for two parameters at the microscopic level implies a significant change in the 
macroscopic limit: the possible loss of hyperbolicity in some regimes. The introduc- 
tion of a non-isotropic kernel of observation, without this dependence on the local 
density, is not sufficient to imply a strong difference of behavior for the continuum 
model. However, it enhances some properties, such as the difference between the 
velocity of the fluid and the velocity of the perturbations of the orientation. 
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It is important to note that the method introduced in [15] works to derive the 
macroscopic model. In particular the concept of generalized collisional invariants is 
still vahd, with some adaptations, and we get the same macroscopic model, except 
for the definition of the coefficients. 

Some questions are left open. The limit here is formal, and we are still looking for 
an appropriate functional framework to obtain more precise results of convergence. 
The rigorous derivation of the mean-field limit of the dynamical system of particles 
is also part of our future work. 

Finally, the next step to this study consists in numerical simulations, in order to 
see how the difference between C2 and Ci can be observed in simulations of the discrete 
dynamical system, or how the particles behave in the regions of non-hyperbolicity. 

A Proof of some statements for section 3 

A.l Expansion of the local density and orientation 

We recall the expressions of Uj'^ and ff: 

u'ix,u,t)= ' (A.l) 

\J^[X, UJ, t)\ 

J%x,uj,t)= f K(^,^-u)vf{y,v,t)^dv, (A.2) 

P%x,u,t)= I t)f|dt;, (A.3) 

Lemma 3. We have the following expansions: 

uj'{x, UJ, t) = n'{x, t) + ea{u- V^) t) + 0(£^) , 
j5%x,uj,t) = p^{x,t) + £au ■ VxP%x,t) + 0(£^) . 

where the constants a and a depend only on the observation kernels K and K, and 

n'{x, t) = ^igli , with f{x, t)= f V fix, V, t) dv , 

p'{x,t) = / r{x,v,t)dv. 

Proof. After change of variable y = x+eC,, let us expand / at first order in e in (A.2). 
We get 

J%x,u,t)= f Km,^-u)v{f%x,v,t)+eC-VJ%x,v,t) + 0{e'))d^dv. 
We have to compute 

Ko{u)=[ K(|e|,|-a;)de and K,{u) = [ Km,^-u)CdC 
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For any rotation R, the change of variable ^ = R{C) gives on one hand 

Ko{uj) = Ko{R{uj)) , 

and so Kq does not depend on u. On the other hand, we get 

i?(KiH) = Ki(i?(u;)), 

which shows that Kiiu) is a vector invariant by any rotation which let uj invariant, 
so it is parallel to u. Given a vector e of §, we have Ki{e) = kie. Then taking one 
rotation mapping u to e, we get R{Ki{u)) = Ki{e) = kie = R{kiu), so finally we 
get Ki{uj) = kiu for all a; G §. Let then a = and we have 

J^(x, CJ, t) 



V {f{x, v,t) + eau- V^fix, v, t)) dv + 0{e^) 



= fix, t)+ea{u- V,.)/(x, t) + 0{e^) . 
Putting this expression into (A.l), we get 



\f{x, t)\' + 2Ea fix, t)-{uj- V,)f{x, t) + 0{e') , 



so 



J^(x, u, t) 



and finally 

u^{x, u, t) 



1 - 



e a 



\f{x,t)\ \ \f{x,t)\ 



■f{x,t)-{u-V,)f{x,t)]+0{e') 



f{x,t) 



+ ea 



But we also have 



\3'{x,t)\ 

{u-V,)f{x,t) f{x,t) {u-V,)f{x,t) f{x,t) 



\3^{x,t)\ \3'{x,t)\ \f{x,t)\ \f{x,t)\^ 



, ^ N^./ X {^-^X^fix^t) ( ^ ( 1 



uj-V.,)f{x,t) 1 



f{x,t) 



{{{u-V.)fix,t))-f{x,t))f{x,t) 



\f{x,t)\ \f{x,t)\^ 
Therefore 

u'{x, UJ, t) = n'{x, t) + ea{u- V^) n''{x, t) + 0{e^) , 

and this is the first part of the lemma. 

After the same change of variable y = x + and expansion in (A. 3), and using 
the same techniques, and the normalization condition (2.5), we get 

p'{x, u,t)= f fix, v,t) + e Ki{u) ■ VJ'ix, v, t) dv + 0{e^) 
= p''{x,t) + eauj- Vxp''{x,t) + 0(5^) . 



This is the second part of the lemma. 



□ 
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A. 2 Proof of Proposition 3.2 

We have to compute (Id — Q ® Q) X, where 

X = [ {{dt + u-\/,){pM,n)+aP{pM,n)+aP{pM,n))K{uJ-n)uduj. 

For convenience, we will write ly, d for i^(p), d{p) in the following. We first give some 
useful formulas to work on the unit sphere. For V a constant vector in M", we have: 

V^{uj-V) = (Id-uj^uj)V, 
■ ((Id -uj®uj)V) = -{n-l)uj-V. 

Then we have that for any constant matrix A 

Vuj ■ ((Id - a; (g) uj)Auj) = A : {Id - nu u), 

where the notation ":" denotes the "contraction" of two operators (if A = (Aj) 
and B = (Bij) then A : B = J2i,j=i,...,n^ijBijj this is the trace of AB^). This can 
be shown when A is of the form Vi ^ V2, using the previous formulas, and then 
extended by linearity. 

We recall the definition of M^n, given in equation (3.5): 

We get, writing cos^ for u ■ Q, and using the notation {■)mk given in (3.7), 

V^M^n = «(Id -u(g) uj)nM^n, 
VqM^u = K(ld - ® Q)ujM^u, 
d^M^n = (cos 6* - {cos 6) mJM^q. 

Using the chain rule, we then get 

{dt + uj- V,)(pM«c) = (1 + (cos^ - {cos9)Mjpk)M^n{dt + u ■ V,)p 

+ pK(Id - (g) n)uM^n -{dt + uj- Vx)^, 

where k is the derivative of k with respect to p. Since Q is of norm 1, we have 
that {dt + u! ■ Vx)^ is orthogonal to Q, and the term f2 ® f2 vanishes. We get 

{dt + u- Vx){pM,n) = (1 + (cos^ - {cos 9) m J pk)M,nidtp + u ■ V,p) 

+ pkM^q{uj ■ dtQ + UJ ^oj : V^^fi), 

where V^f^ is the gradient tensor of Q that is to say {'Vx^)ij = dx^^j- We then have 

P{pM^n) = i^{p)V^ ■ ((Id -uj0uj){{u- Vx) n)pM^n), 

= pu{p)[nn ■ (Id - w ® w)((V^fi)^w) + Vo; ■ ((Id - a; ® uj){\/x^)^uj)]M^a, 
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where the notation denotes the transpose of operators. Hence, using the fact 
that {'Vx^)'^oj = (u-Vx) ^ is orthogonal to Q, and the formula given in the beginning 
of this section, with A = (Vx^)'^, and we get 

P{pM^u) = pu{p)[-Kcos9uJ^u : {Vx^f + {Vx^f : (Id - nu ^ uj)]M^n, 
= pi^[Vx ■ — (n + KCOs6) u ®uj : V^f^jM^Q. 

Similarly, for the operator P, we get 

P(pM«n) = />(p)V^ ■ ((cu ■ VxP) (Id - a; ® uj)VtpM^^) 

- d{p)V^ ■ (ipM«f^(Id -uj®u)Vxp+{oo- VxP) V^pM^n), 
= p{p - Kd)V^ ■ ((Id - a; ® ® VxP)coM^n) 

- lpd[Vu^ ■ ((Id - a; ® uj)Vxp) + kQ ■ {Id - u ^ u)VxP]M^q. 

But we have u = nd, so i) — nd = dk. And we have 

■ ((Id -u0u){n® VxppM^n) 

= Vuj ■ {{Id - cu IS) uj){^l VxP)co)M^n + ■ {Id - u uj){n ® Vxp)^^] 
= [{n (g) Vx'P) : (Id -nuj®uj) + k{1 - cos^ 6) to ■ Vxp]M^n 
= [Q ■ VxP + {i^ sin^ 9 — n cos 6) u ■ VxPIM^q. 

Hence 

P{pM^n) = pdk[Q ■ VxP + (/« sin^ 6 -ncos6)uj ■ VxP]M^n 

+ ^pd[{Kcos9 + n — l)u ■ VxP — ■ VxPjMj^^. 

Finally we can write X = Xi + X2 + X3, where 

Xi= h^,{cos6)'yi{cos9) u M^,ndu, 
^2 = h^{cos6) u ® cj(72(cos 9) VxP + pi^dtQ) M^ndu, 
= hf,{cos9)ys{cos9) u{uj S u : Vx^) M^ndw, 

with (using the notation ci = (cos^)a/«) 

7i(cos6') = (1 + (cos 6* — ci)pk)dtp + apvV x ■ ^ + ap{dk — ^dK)Q ■ VxP, 
72(cos^) = 1 + (cos^ — ci)pk + ap{dk{Ksm'^ 9 — ncos9) + \dK{Kcos9 + n — 1)), 
73 (cos 6') = pK — apv{n + kcosO). 

To do the computation we write u = cos9Q + sin9v, with v G S„_2 (identified 
with the set of unit vectors which are orthogonal to Q). We take the following 
convention: /„g§^_2 d^; = 1, and we have 

/ a{u)du =^ r f a{9, v) sin"-^ 9 dv d9, 

I vdv = Q, and j v®vdv = — — (Id — 1] (S) f2) , 
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where Vn is a normahzation constant (we will not need it in the following). These 
results are still vahd when n = 2, with J^^^^ dv = ^{a{vo) + a(— wq)), vq being one of 
the two unit vectors orthogonal to Q. Using these formulas, we get 

7(cos^) M^f^wdo; = (cos^7(cos^))m^ il, 



u 0uj-f {cos 9) M^ndu = (cos=^ ^7)^/,^ ® + -^^^^^^-^^^(Id - ® fi). 
/wgs n — 1 

So we have (knowing that dt^l is orthogonal to Q): 

{id-n® n) Xi = 0, 

n — \ n — 1 

To compute (Id — Q X^, we first remark that 

(Id - (g) Q) uj{u (g) w : V^fi) = sin 9v{u ■ {u ■ Vx)^) = sin^ 9v{v ■ {u ■ Vx)^), 

since (u-Vx)^ is orthogonal to Q. But we have /t,g§^_2 v{v^v:'Vx^) dv = 0, because 
the integrand is odd with respect to v, and then we get 

(Id - 1] ® X3 = (sin^ 9 cos 9 73 h^) v0vdv{Q- V^)fi 

{sm^9 cos 6* 73 /i«)a/. 



n — 1 

since {Q ■ Vx)^ is orthogonal to Q. 

So we have that (Id — ® f2)X = is equivalent to 

pK{si\)? e h^) M^dtVt + (sin^ cosO -^^h^) ' x)^ + (sin^ 6'72/Ik)m«Vx-P = 0. 

For any function 7(cos 9), we denote by (7(cos 9)) the mean of 7(cos 9) follow- 
ing the "weight" sin"^ 9 h^{cos9)Mfj^Q, that is to say 

/(77(cos^)/i^(cos^)e^'=°^^ sin"^d^ 

(7 (cos ))^^^ - J- h,(^cos9)e-^°^<^ sin"^d^ ' 

We have 



(sin^ ^ ^k(cos ^)7(cos 6')) 



and so, dividing by ^{5111^ 9h ^{cos 9)) we finally get that (Id — f2 ® Q)X = is 
equivalent to 

p {dtn + C2{n ■ Vx)^) + A (Id - ® n)\/xp = o, 

where the coefficients are given by 

C2(P) = ;^(cos^73(cos^))^~^_^ and A(p) = ^{-f2{cos9)) 
We finally get, writing ci = (cos 6*)^ , 

C2 = Ci — a d {nci + K (cos^ ) ; 

A = ■^+p^[(ci — Ci + ad («;(sin^ 6')^-;^ — ri Ci) ] + | a d (kci + n — 1) , 

which are exactly the expressions given in equations (3.19)-(3.20), and this ends the 
proof of Proposition 3.2. 
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B Asymptotics of the coefficients 
B.l Proof of Proposition 5.1 

We recall that the two linear operators L and D on the space of polynomials are 
defined by 

L{P) = -(1 - X^)P" + (n + 1)XP' + (n - 1)P 
D{P) = -{1-X^)P' + XP. 

We first give a preliminary lemma which will be helpful to construct the polynomi- 
als Hp and . 

Lemma 4. Definition of the polynomials. 
Let Q be a polynomial and G N. Then 

• There exists one unique polynomial P such that L{P) = Q. 

• There exists one unique polynomial of degree at most N such that 

D{P^){cos9) = Q{cos9) + 0(6'2(^+^)) as 9 ^ 6. 

Proof. For the first point, if the leading term in a polynomial P is akX^, with ^ 0, 
then the leading term in L{P) is [k{k — 1) + k{n + 1) + {n — l)]akX'' , and so L{P) ^ 0. 
So the linear operator L is injective from Mp[X] to Rp[X], and therefore it is bijective. 
For the second point, the idea is to remark that 



- Xf) = {2k + 1)(1 - X)'' + {k + 1)(1 - Xf+\ 

so we write the polynomials in the basis {(1 — X)'', k G N}. We get that a poly- 
nomial R is such that R{cos9) = 0(9'^^^'^^^) if and only if, in this basis, its first 
coefficients up to order {X — 1)^ are zero (because 1 — cos 6* = ^9"^ + 0(6*^) in the 
neighborhood of 0). We write Q = Er=o^fc(l - ^)'' and P^ = Ef=o«fc(l - 
and we get that 

D(P^)(cos^) = Q(cos^) + 0(^W)) ^ h = bo , ^, , 

[{2k + l)ak - kak^i = bk VA;g|1,7V|. 

Since this induction relation defines in an unique way the coefficients ak for k G 
|1, A^], this ends the proof. □ 

With this lemma, we can now define the following sequences of polynomials Hp 
and Gp , this last ones being of degree at most N — p: 

'l{Ho) = 1 fD(G'f ) (cos ^) = 1 + 0(^2iV) 

L{Hp^,) = -D{Hp) \{D{G^^,) + L(G,^))(cos^) = 0{9'(^~^)). 

Since the operator L is odd and D is even, it is easy to show that the polynomials Hp 
have the same parity as p. If we express the operator L in the basis {(1— X)^, k G N}, 
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we are able to get the induction relation for the coefficients of the polynomials in 
this basis. We have 

L((l - X)'') = (n + - 1)(A; + 1)(1 - X)'= - k{n + 2k-l){l- X)''~\ 

So we have 

p N—p 

Hp = Y. - X)\ and G^, = ^ ^(1 - X)^ (B.l) 

fc=0 k=0 

where and are given by the following induction relations for (with the conven- 
tion that t^_^_i = Jf^i = a^i = 0): 

b'o = ^,, 

(n + k- -in + 2k + = '-tillf, - ^Ci, Vp G N, Vfc = [0,p + Ij. 

k (2/c+l)(2A:-l)...3 (2fc+l)! ' 

fS^af' - i^,at\ = {n + k- l)al - {n + 2k + l)af+„ G N, Vfc G N. 
We define then the remainders R^q and R^^o by 

N N 

<o(/^) = Kill) - E Hpili)^', <oo(/^) = hM - E 

p=0 p=l 

It is an easy matter to see that, for a given polynomial P, we have 

sin^<9e(P(cose)sine) = D(P)(cos^) sin^, (B.2) 
-9e(sin"-2 6'ae(P(cos^) sin^)) + (n - 2) sin"-^ ^F(cos^) = sin"~^ ^L(P)(cos 

and then we get 

Z*(P(cos^) sin^) = {L{P) + kD{P)) (cos 9) sin 9, 
where the operator L* is defined in (3.10) by 

L:gi9) = -sm'--ee-^^°^'^,ism--'9e^^°^'g'{e)) + ^,g{9). 
Since we have by definition L* (/ik(cos 6^) sin 6') = sin 6*, we get 

N 

L;(i?^o(cos^) sin^) = sin ^ - E(^(^p) + «^^(^p))(cos^) sin^K^ 

p=0 

= -k^+^D{Hn){cos9) sin0, 

N 

^:«oo(cos^) sin^) = sin^ - Y.iL{G^) + kD(G'^^))(cos ^) sin^K'^ 

p=i 

N-l 

= -L{G%){cos9)sm9K-^ + E ^^-"'0(9^^''-^^) sm9. 

p=0 
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To get estimations for the averages of the form (/(6')i?^^(cos 6') sin^ 6')Af« (with e 
standing for or oo) we first remark that, for a function g belonging to the space V 
(a "weighted Hi'') defined in (3.12) by 

V = {g\{n- 2)(sin^)t-2^ g ^^(0, n), {sm9)^~'g G H^iO, n)}, 

we have the following Poincare inequality: 

{g{9)L:g{e))M. = {g\ef)M. + (n - 2){^^g{9f)M. ^ {n - 2)((^(^))^)m.. 

Hence, for n ^ 3 and g & V , using Cauchy-Schwarz inequality, we get that 

{g^M. ^ ^^{9{ef)MA{hg{9)Y)M.. 

Since g^iO) = hi^{cos9) sin 6* belongs to V, we get that g^i;{9) = R^^{cos9) sin 9 also 
belongs to V. 

We are now ready to do the estimations. For / such that 9 t— ?■ f{9) sin a' 9 belongs 
to L^(0,7r), we get, using Cauchy-Schwarz inequality, 

\{f{9)R^^,{cos9)sm'9)Mj < ^ {{R^^oi^os 9))^ sin' ^)a/.(/(^)^ sin^ 9),,^ 
^ ^ V'((L*<o(cos 9)y sin' 9)Mj{fi9y sin' 9)m. 

^ ^K^^'^mHj,)icos9))' sin' 9)mJ if (9)' sin' 9)m.. 

Hence using the expansion as k — > given in (5.1), we get the final estimation: 

(/(6')i?^o(cos 9) sin' 9) = 0{k^+^) as k ^ 0. (B.3) 
Similarly, if |/(^)| = 0{9'^) in the neighborhood of 0, using Lemma 2, we get 

\{fmloo{^os9)sin'9)Mf ^ ^^(/(^)'sin2^)M. 

X /[L(G'j^)(cos 9)k-'' + «;""0(^2(^-P))]2 sin^ ^\ 
^0{k-'^-^)xO{k-"'-^), 

which gives 

{f{9)R^^^{cos9) sin' 9)m^ = 0{k-^-^-^) as k ^ oo. 

Now, since we have the expression (B.l) of the polynomials Gp , we get, by 
definition of -R^oo; 

p=0 
N+l 

p=0 
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Since (1 — cos 6*) = 0(6 ), we finally get, using Lemma 2, 

N+l 

(/(^Xoo(cos^)sin2^)M. = {mR^X'(^ose)sm'e)M. + E 

p=0 

This ends the proof of Proposition 5.1, in the case n ^ 3. 

We suppose now that n = 2. The case k — )• is easy, since we have the following 
Poincare inequality: 



We get the same estimations, replacing (ri — 2) by e 



-2k. 



I (/(eXo(cos ^) sin^ 9)^,^ I ^ e2''«:^+ V((^(^^) (^^^ ^))' ^^^^ ^) (/(^)' ^^^^ ^) ' 

which gives the estimate (B.3) since e^'^ = 0(1) when k — )• 0. 

The case k — oo is different, since we are not able to get a better Poincare con- 
stant. But we have an explicit expression for = h{cos9) sin 6', given by (3.13): 

£ _ 7r /o^--^°-^dv9 

It is also easy to see that, in this case, the coefficients appearing in the defini- 
tion (B.l) of the polynomials , are zero when p ^ 1. Therefore we get that = 
for p ^ 1. 

We have ^(G'f )(cos^) = 1 + 0{9^^), so with the formula (B.2), we obtain 

^^(^^(cos^) sin^) = l + 0{9^^), 

so we get that G^ {cos 9) sin 9 = 9 + 0(9'^^'^^) since 9 Gi{cos9) is continuous 
as 6' — 0. Actually, this is the Euler formula for arctan: if we write t = tan | 
we get 1 — cos^ = -^s, sm9 = and then, using the expression (B.l) of the 
polynomials G^ with au = '^2k+iy. ' obtain 



f N r,2k(y\2 f2k 

Now, using the explicit expression of gt-, we have 

i?^^(cos^) sin^ = ^,(^) - Gf (cos^) sin^«;"^ 

= ^-1(0 _ Gf (cos^) sine) - - {{"'^'"'^^^ 
= K^^O{9^''+^)-r^{9). 



We have 



^2 g— KcosS 
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and so, using the estimate (5.3) with n = 2, we get 



4 fTT -2re cos 6»K cos 6» J /Q 4 

\l K I /M. ^ ^2 (J-g_,,os^d^)2 J-g.eose^J^ - ^2 ,os ed^)2 ^ ^ 

Therefore, using Cauchy-Schwarz inequahty and Lemma 2, 



|(/(e)r-(^)sine),,J ^ v'(/(^)'«i^'^)M«((^^(^))')A/. = 0{K-^~'e~'n 
so we get the final estimate 

(/(^)^^oo(cos^)sin2 6')A/, = 0{k-^-^-^) as -> oo, 
and this ends the proof of Proposition 5.1. 



B.2 Tips for the general case 

Here we give some tips to perform an asymptotic study of the coefficients when u 
depends also on u ■ Q. 

We will have to take averages against functions of the form 9 gK(p,cose)^ where 

1 r^* 

We want to get for example an expansion as the noise d is large or small. So we are 
only interested in the dependence on cos^, and we will drop the dependence on p 
for clarity. We suppose that the function 6 i— )■ z/(cos 9) is positive, smooth, bounded 
below and above, and we introduce the parameter k = ^, trying to expand with 
respect to k. We write = jQu{x)dx, so we have = Ka{ji). 
The first step consists in the expansion of {f{0))^j given by 

^ /o"/(g)e-(-^)sin"-^gdg 
Wl ))m^ /(7e''-(=°^«)sin"-2 0d^ ' 

As before, we can easily do a Taylor expansion when k — )■ 0, and we get a result 
similar to (5.2) involving quantities of the form /J^ f{9)a{cos6yd9. Unless we know 
explicitly cr, we cannot say anything interesting. 

When K ^ oo, the strategy is the same: we do the change of variable, setting t = 
cr(l) — o"(cos^)), and a{t) = a~^{a{l) — t), where is the inverse function of cr 
(which is increasing since u > 0, actually we have a{t) = cos 9). We get: 

^Ka{cose)-n~2n.n_^Ka{i) ^ /(arccos a(t))e- ^* 



j.mg«<x(cosf^)g.^n-2^^^ _ gK<x(i) / _iy y n 

^ Jo u{a{t)){l - a{ty)^ 

where T = cr(l) — a{—l). Since a(0) = 1, if we know the expansion of / around 
and u in the neighborhood of 1, it only remains to get a Taylor expansion of a 
around to use Lemma 1 (Watson's Lemma). We can compute the derivatives of a 
by induction. We have 

a'it) = ^ 



u{a{t)) 
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and this gives immediately a^"-* = F„(a(t)), with the following induction relation 
for F„ 

This gives us the Taylor expansion of a at t = up to order A^: 

aW = l + E/n(l)3 + 0(t^^')- 

n=l 

Since we have a'(0) < 0, it is possible to get the analogous of Lemma 2: for a 
function / such that f{e) = 0(9^^), then {f{0))-^ = 0{k~^) as k ^ oo. Finally, 
we get an expansion of {f{0)) q which depend only on the first derivatives of z/ at 1 
and on the local behavior of / around 0. 

The second step consists in expanding {f{6)gi^{6) sin 9) £^ , where L*Jj^{6) = sin 9, 

the operator L* being defined by 

L:g{9) = -sin2-"^e-'^'^(^°^^)A(sin-2^e«.(co..)^,^^)) ^ i^AO)- 

It is easy to see that we have 

Ll{P{cos9) sm9) = {L{P) + {cos 9) D{P)) {cos 9) sm9, 

so if we set ^k(^) = hf^{cos9) sin 9, we can decompose in a way similar to Propo- 
sition 5.1: 

TV 

K{cos9) = Hp{cos9)kP + R^^oicos9), 

p=0 
N 

K{cos9) = Y.G^{cos9)k-p + R^^^{cos9), 
p=i 

where Hp are the functions (not necessarily polynomials) and the polynomials 
of degree at most N — p given by the following induction relations: 

'l{Ho) = 1 

L{Hp+M = -u{f,)D{Hp){i,) 

(D{Gf){cos9) = 1 + 0(^2JV) 

[iD{G^^^) + ^)L{G^mcos9) = 0{9'(^-^)). 

Again, the remainders satisfy the following estimations, for any function / such 
that 9 ^ f{9)sm^9 belongs to ^^(O,^) and such that |/(^)| = 0{9^(^) in the 
neighborhood of 0: 

{f{9)R^Q{cos9)sm'^9)Q^ = as k ^ 0, 

{fi^)Kooicos9)sm^9)-^ = 0{k-^-^-^) as k ^ oo, 

and this allows to get an expansion of {f{9)hK{cos9)sm'^9)g when k — j- and 
when K — ?■ oo. 
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